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We calculate the collective mode spectrum in models of superconductors with attractive interac¬ 
tions in s-and d channels as a function of their relative strength, across a phase diagram that includes 
transition between s, s + id, and d- wave ground states. For one-band systems, we recover the largely 
known results for the phase, amplitude, and Bardasis-Schrieffer modes of pure s or d-states, and 
show how the well-defined Bardasis-Schrieffer mode softens near the s + id phase boundary and 
evolves in a characteristic manner through the s + id phase as a mixed symmetry mode. For two- 
band systems, we consider a model of hole-doped Fe-based superconductors, and find in the case 
of an s-wave ground state a well-defined Bardasis-Schrieffer mode below the lowest gap edge, as 
well as a second, damped mode of this type between the two gap energies. Both modes soften 
as the s+id phase is approached, and only a single “mixed-symmetry Bardasis-Schrieffer mode” 
below the pairbreaking continuum propagates in the s + id phase itself. These modes coexist with 
a damped Leggett mode with collective frequency between the two gap scales. In the pure d-state, 
no Bardasis-Schrieffer type s excitonic mode exists at low T. We briefly discuss Raman scattering 
experiments and how they can be used to identify an s + id state, and to track the evolution of 
competing s- and d- interactions in these systems. 


I. INTRODUCTION 

In contrast to cuprate superconductors, which are be¬ 
lieved to have universal d -wave symmetry*^ the Fe-based 
superconductors (FeSC) appear to manifest generic s- 
wave pairing, with an order parameter that likely changes 
sign between Fermi surface (FS) sheets^— On the other 
hand, calculations of pairing by exchange of spin fluctu¬ 
ations have suggested from the early days of research on 
these materials that the d -wave channel can be strongly 
competitive, and might under some circumstances be¬ 
come the dominant pair symmetry^i Within an RPA 
treatment, it was argued that overdoping either by holes 
or electrons away from a 6 electrons/Fe parent mate¬ 
rial should lead to a d -wave ground state<&£ This gave 
rise to the possibility that the symmetry broken s-wave 
phase could make a low-temperature transition to a d- 
wave phase, and several authors argued that this should 
proceed through an intervening s + id state, in which s 
and d wave functions are combined with fixed relative 
phase 7t/2 and time reversal symmetry T is brokem 10 : 11 
While such a transition would be of great potential inter¬ 
est, and represent the first example of its kind, an s + id 
state is not trivial to detect. Although T is broken, the 
state is not chiral, and thus does not manifest sponta¬ 
neous edge currents as discussed, e.g. in the context of 
the p + ip' state of Sr 2 Ru 04 . Its quasiparticle excitations 
are fully gapped, but so are those of the s wave state out 
of which it evolves, so thermodynamic signatures of the 
transition are likely to be weak. 

The proximity of different pairing channels is an un¬ 
usual situation in superconductors(SC), but was stud¬ 
ied quite early in the pioneering work of Bardasis and 
Schrieffer,— where the effect of fluctuations in a subdomi¬ 
nant pairing channel was investigated in a conventional s- 
symmetry ground state. The motion of the order param¬ 


eter was found to include a collective mode correspond¬ 
ing to the oscillation of the phase of the subdominant 
pairing channel, with the q —> 0 frequency dependent 
on the difference between the inverse of the two pair¬ 
ing interaction components. This frequency is located 
below the pair-breaking edge of the condensed s-wave 
system. This “Bardasis-Schrieffer (BS)” mode (some¬ 
times referred to as a particle-particle exciton) was never 
convincingly observed in conventional superconductors 
owing, presumably, to a dearth of systems exhibiting a 
strong s-d competition. Taking a hint from the proximity 
of s and d channels predicted by spin fluctuation theory, 
FeSC can be viewed as excellent candidates to probe such 
modes. A search using Raman scattering was proposed 
in FeSC by Devereaux and ScalapinoJ^ Recently, two 
electronic Raman measurements on Bai- a ,K a; Fe 2 As 2 14115 
and NaFei-aCo^As^ found features associated with BS 
modes. The exact identification of these features with a 
BS mode is hindered by the fact that these systems posses 
multiple gaps and the exact nature of possible collective 
modes and their evolution across a typical doping phase 
diagram is not clearly known. Nevertheless, these dis¬ 
coveries raise the prospect of systematic studies of the 
interaction strengths and collective modes in different 
channels in FeSC for the first time. 

To facilitate the interpretation of experimental results 
in this area as they develop, we provide here the col¬ 
lective mode spectrum in a simplified model of an un¬ 
conventional superconductor with competing pairing in 
s and d channels. The prime focus of this work is the 
study of collective modes in a non-s wave ground state. 
We thus explicitly account for the possibility of a tran¬ 
sition to an s + id state between pure s and d phases 
and investigate collective modes in the s + id state and 
report a mixed symmetry collective mode which couples 
the amplitude and phase sectors of the fluctuations in 
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the SC order parameter and exhibits oscillations in both 
s and d symmetry channels. As a demonstration of our 
approach, we also study the simpler, one-band case, and 
drive it through s to s+id to d transitions. We reproduce 
the well known results for collective modes in single band 
superconductor and use them as a benchmark to discuss 
the differences in a multiband system. We demonstrate 
the existence of the mixed symmetry collective mode in 
both one and multi-band cases. 

The study of collective modes in ‘exotic’ non-pure- 
s wave superconductors has some history: Hirschfeld 
et alii studied the analog of the 3 He-B “squashing 
modes” in the Balian-Werthamer p -wave ground state, 
discussing how such modes could be observed in opti¬ 
cal conductivity. Wu and Griffinii studied the possi¬ 
bility of an s-excition for systems with d -wave ground 
states, and showed that this mode does not propagate 
at low T. More recently, in the context of the FeSC, 
Devereaux and Scalapino investigated the role of a BS 
particle-particle exciton on the Raman spectrum of an 
s± superconductor^ Khodas et. al. studied Raman sig¬ 
nature of this mode and the role of density fluctuations 
in Fe-selenidesJ^ Bittner et. al. derived general expres¬ 
sions for the collective modes in the s-wave ground state 
for noncentrosymmetric systems^ 

Works on collective modes in T-broken superconduc¬ 
tors are somewhat rarer. Several papers generalized the 
well-known collective modes of the 3 He-A phased! to 
p + ip' and d + id' superconductors in situations where 
the two harmonics in question corresponded to the ba¬ 
sis functions of a 2D representation of the symmetry 
group ; 22 ’ 23 Balatsky et. al* A phenomenologically dis¬ 
cussed a ‘clapping’ type orbital mode analogous to those 
discussed in Ref. [ 23 , but for a general situation where the 
two harmonics d and d' were not necessarily degenerate. 
The same mode and its detection in Raman experiments 
was also discussed in Lee et. ali^ 

This mode corresponded to an oscillation of the rela¬ 
tive phase of the two components d x 2_ y 2 (d) and d xy (d'), 
two distinct representations of the tetragonal group. 
More recently, one of the authors^ and Marciani et. 
studied, in the context of FeSC, the Leggett modes (the 
oscillations of the relative phase between order param¬ 
eters on two bands: see Ref|28j) in a special s-wave T 
broken s + is' SC where it was found that the mode soft¬ 
ened at the boundaries of the s + is' state (see also Refs. 

iffl. 

The current work is more along the spirit of the last 
work above in the sense that the collective modes in a 
multiband s + id SC are investigated. The mixed sym¬ 
metry collective mode that we find also softens at the 
boundaries of the s + id state. Besides the mixed sym¬ 
metry collective mode in the s + id ground state, we 
also find, in the s —wave ground state, a damped Leggett 
mode (which has s-wave symmetry) residing between the 
multiple gaps in the system (in our language we do not 
consider them to be true collective modes of the system, 
but nevertheless appropriate response functions will show 


broad peaks). Additionally, besides the usual BS mode 
residing below the minimum gap in the system, we report 
another damped BS mode residing between the multiple 
gaps. Within our model, we can show that there is only 
one BS mode below the minimum gap for any interaction. 
The presence of multiple BS modes is therefore expected 
to be a generic feature of multiband systems. 

In this work we have ignored the coupling to the den¬ 
sity fluctuations. The main qualitative effect of ignor¬ 
ing the density fluctuations is related to the Bogoliubov- 
Anderson-Goldstone(BAG) mode. Although the BAG 
mode is expected due to spontaneous breaking of U(l) 
symmetry during the SC transition, coupling to charge 
density fluctuations implies that the oscillation of the 
gauge degree of freedom is identical to the usual plas- 
mon (“Anderson-Higgs mechanism”) ^ In conventional 
cases of 1-band and 2-band SC, it is well known that 
coupling to density fluctuations does not affect the mass 
of the Leggett^— or the BS modest Since the argu¬ 
ments are based on gauge invariance and symmetries, we 
assume without proof that same hold for our multiband 
system. We expect our work to provide useful insight in 
terms of number of collective modes to be expected in a 
system and detecting a non-trivial multiband SC ground 
states (s + id is our case). Although there are many ways 
to model the multiband scenario, we limit our consider¬ 
ations to a minimal model that can be readily applied to 
FeSC, as will be described later. 

We study the collective modes by studying the pos¬ 
sible excitations in the system within linear response in 
different angular momentum channels. This method is 
sometimes referred to as a generalized random phase ap¬ 
proximation, and is known to yield results identical to 
those obtained from the kinetic equation method^!. We 
explicitly derive a 1-band case and extend the formula to 
the multiband scenario. We stress that the formulation 
has the great advantage of spitting out all the collective 
modes in a clean SC in all angular momentum channels 
with minimal effort. The biggest advantage is its scalabil¬ 
ity to multiband or multiorbital systems. The rest of the 
paper is organized as follows: In Sec. El we specify our 
1 band model, derive the collective mode equation and 
study the collective modes; reproducing the well known 
results. In Sec. EH we discuss the modeling of a FeSC 
with 3 pockets, discuss the collective modes and highlight 
the differences with 1-band model. In Sec. IIVI we discuss 
our results in connection to FeSC and some of the recent 
Raman experiments. We summarize our main findings 
in Sec. The Appendix presents details of some of the 
calculations. 


II. COLLECTIVE MODES IN A 1-BAND 
MODEL 

We revisit this simple model as this helps us in two 
ways: we can systematically tune the system through a 
s to s + id to d transitions and trace the collective modes 
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across the phase diagram (through the s+id phase whose 
collective modes have not been addressed before); we will 
then use this result as the benchmark against which the 
multiband case will be compared. Since we are interested 
in studying the collective modes in the SC state, we only 
retain the interactions in the particle-particle (p-p) chan¬ 
nel. As discussed in the introduction, the interactions 
in the particle-hole (p-h) channel (which couples the SC 
fluctuations to density fluctuations) will be dropped. The 
logic of our presentation will be the following: we start 
with an s —wave ground state; use the d —wave interac¬ 
tion as our tuning parameter to generate a phase diagram 
that scans through the s, s + id : d wave regions; we then 
find the collective modes with both s and d symmetries 
in each region and consistently track them as the d —wave 
interaction is tuned. 


A. Model and phase diagram 


Our one band model is a 2D Fermi liquid (FL) with 
the interaction V(p, —p;k,—k) = V(p,k) in the pairing 
channel. We decompose this interaction into different 
singlet angular momentum channels (limiting ourselves 
up to the d —wave harmonic): 

V(p,k) = U s + U d f u fp, (1) 


where /g = \/2cos2Op. All the vectors are by definition 
on the circular fermi surface (FS) and Op is the angle of 
k measured from the k x - axis. We follow the convention 
where repulsion is denoted by the positive sign of E/’s. 
We keep U s fixed and tune U d . Within weak coupling, 
this results in the usual self consistency relation for the 
order parameter Ap££ 


A p = 



( 2 ) 


where f K stands for TJf n f pj) 3 • Taking as input that 
the only stable solutions are s, s + id , and d states (no 
s + d)r^ and writing Ap = A s + A d fp (with A s ’ d as 
constants) we arrive at 


A s 

A d 



A s 

u*+el+\Atf’ 
ul+el + \A % r 


( 3 ) 


Since we work in the FL regime, we implement ft = 
1 / 2 d f If-cfe which leads to the definition of two dimen¬ 
sionless parameters u s = V 2 dU s and Ud = V 2 DU d where 
V 2 d is the 2D density of states at the fermi surface. This 
model is then easily solved (see Appendix A) and the 
resulting phase diagram is schematically plotted in Fig. 
H|i. 


The quantities that change with Ud are T c , A s , and 
A d . To remove the energy cut-off (A) dependence of our 
results, these quantities are normalized to Ag, the gap 
value at Ud = 0 (the pure s —wave state). We thus work 
with the normalized parameters: 
a s = A s /Ag, ad = A d /Ag, 77 = ad/a s . 

The most relevant points that define the boundaries of 
the s + id phase are at T = T c and T = 0. T c (obtained 
by setting q s j 4 0 ) across the phase diagram is given 

by 


In 


2 7 A 

nT r 


1 


1 


= min <-,-> > 0 , 

u s u d 


( 4 ) 


where 7 ~ 1.78. The boundaries of the s + id state can 
be found after rewriting the self consistency equations as 
(see Appendix A for details) 



and setting 77 — 0 and rj —> 00 in the second equation. 
This results in the s/s + id boundary at u s = Ud and the 
s + id/d boundary at Ud = 2u s /{2 + u s )^ Fig. [TJd shows 
the calculated gaps a St d as Ud is tuned through the s + id 
state. 


B. Collective modes-formulation and results 

We employ standard linear response to study the col¬ 
lective modes in this system. We provide a simple deriva¬ 
tion for the 1-band model as this will help us generalize 
the formula to the multi-band case with ease. The Hamil¬ 
tonian in the SC state is given by 

H = where ^ = (4 t ’%)’ 

k 

- a s 'n + 

+ = -I>69X c V«>' (6) 

q 

where do is a 2 x 2 identity matrix, <7 is a vector of Pauli 
matrices and R, I stand for real and imaginary part of 
the order parameter. The perturbing fields (originating 
from fluctuations of density and the order parameter) 
that couple to this Hamiltonian have the form: 

SH^t) = (SD^-SA^+SA^v 2 )e—*.(7) 

It is convenient to deal with perturbing fields indepen¬ 
dent of the internal variable k in order to formulate the 
linear response problem. This is achieved by writing 

SAf t = YSAf L ff, 

q,k 9 J k ’ 

L 

s+i = E iD Pf> 

L 


( 8 ) 
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FIG. 1: (a) The schematics of the phase diagram of the 1-band model, (b) The calculated evolution of the s— and d— wave 

gaps a s and ad (normalized to the pure s— wave gap) as a function of the d— wave interaction strength Ud- (c) The evolution 
of the collective modes(solid lines) in a 1-band model as \ud\ is increased. The BS mode softens as the s + id boundary(shown 
in light blue dots) and acquires a mixed symmetry (MS) character in the s + id state and softens again at the other boundary. 
The dashed black line denotes the minimum gap in the system. Collective modes are well defined only below the minimum 
gap. There are no other damped resonances in this model. Here u s = —0.5. 


where L are the different orthogonal angular momentum 
channels. (These will correspond to the irreducible repre¬ 
sentations of the point symmetry group in the presence 
of a lattice). In the usual case of density fluctuations, 
there is no dependence of on k (as it corresponds 

to non-equilibrium fluctuations of total density). This is 
why the higher angular momentum channels, in this sim¬ 
ple model, are not affected by the Coulomb force. This 
will be utilized later. For now, proceeding in complete 
generality, we are then led to 


we can make use of the following relations 
'ido-ni' = c|c| — C-J-CJ,, 
4' t cr 2 4' = -i (cjcj + c t cj.^ , 
'F t cr 3 'l' = cjc t + cjcj., 
to write down, 

k 


= E fZH SA f ,L °i + S K ,L ^) 

k,L 

= (?>-“*. (9) 

L,i 


-2SAl = ~^2 v (q,k)6(^la 2 \tg), 

k 

SD ? = 


(15) 


(16) 


where 

^(q) = E/K-^-’ ( 10 ) 

k 

and the perturbing field 

SF t L = (- SA R ’ L , SA i ’ l ,SD l ). (11) 

The self-consistency equation in Eq. © is then written 
as: 

(A’) 1 = (12) 

k 

where V LL is defined through 

V(k,q) = ( 13 ) 

L,L' 

Starting from 

jaj = 


where V^= 2n e 2 j q. Following the argument around Eq. 
m and abbreviating V LL as V L , we can write Eq. m 
as: 

25Af' L = - V l 6(0f), 

-2SAl’ L = -V L S(0%), 

SDi = VlAipli). 

(17) 

We will work in the limit of q —> 0 so that q will only be 
retained in Dq to account for the singular nature of the 
Coulomb interaction. 

The statement of linear response is that the change in 
an operator due to the applied perturbation is given by 
(recalling that Q = — (’k’l'f)) 

6(n^)(Q) = J2^fj L 'm F J L \Q), 

j,L ' 

n E (Q) = f f%f% Tr [G(K)a z G(K + Q)aj]. (18) 


(14) 
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Combining Eqs. (El)-©, we arrive at the one-band 
system of equations: 

53{n«'' -2[V L ]- 1 6 LL 'S ij }6Ff = 0. 

3,l> L 

(19) 


See Appendix B for explicit form of the mode equation. 
The non-trivial solutions of this set of equations are the 
collective modes of the system. 


1. General considerations 

We see from Eq. m , that the collective modes can 
be found once the interactions V L and the polarization 
bubbles H^ L are known. These depend on the specifics 
of a microscopic model and can be easily computed. In a 
particle-hole symmetric system, quite generally, we will 
have: 


n = 

IM 

~e 2 + A 2 R - A 2 } , 

n£f(^„) = 

!M^- 

^£ 2 -A 2 r + A 2 ] , 

H33 (iQ n ) = 

IM 

+ e 2 -A 2 R -A 2 ] , 

nfg 1, (ifl n ) = 

L v dX [n " A ' ] 

-n 

n 23'(i^ n ) = 

IM™* 

} = 

n (iVL n ) = 

1MM’ ■ 



( 20 ) 

D± = (w m ± ^) 2 + e 2 + |A| 2 . A R J are the ground 
state properties and are taken as input from the analysis 
in the previous section. Although the temperature evo¬ 
lution can be tracked, we shall perform calculations at 
T = 0 as the calculations are tractable and already very 
informative. We make use of the following integrals: 




D+D_ 


IK ly + 


1 




Ik ly + 


- [- 

A° h4E' 


(21) 

1 / 1 

1 

(22) 

(AO ) 3 h4E E2 

V2A °) 


where E 2 = + a 2 + OL 2 d fp and 

r , 2DA°/o(in n ) = / 

Jk 

= [ 

Jk 


V2DA®l4(in n ) = / — 
Jk E 


1 

1 

E E 2 + 

(M' 

V 2A 0 ) 

1 

n 

E E 2 + 

(M 2 

1 

n 


E 2 +{£$) 


2 ’ 


(23) 

(24) 

(25) 


where fg =1; fg = \/2cos20. It is worth noting that in 
the s— wave ground state (ay = 0), I 2 = Iq and I 4 = |/o- 
Analytic continuation to real frequencies is performed by 
ifl n —> tt+i5. As one tunes Ud, the ground state changes, 
and the ’s need to be calculated at every Ud- 

In what follows, we will ignore the coupling of the col¬ 
lective modes to the charge sector as justified in the Intro¬ 
duction. We refer the reader to discussion in Ref. [ 2 ^ and 
references therein where explicit coupling to the charge 
sector, within the same formalism, is presented. The re¬ 
sults from now on therefore pertain, strictly speaking, to 
the collective modes in a ‘neutral’ SC. This simplification 
allows the formulation of the whole problem in a 4 x 4 
space of SAf'J. In the chosen gauge (where the s— wave 
condensate is chosen to be real), R maps onto the ampli¬ 
tude sector and I maps onto the phase sector. We now 
investigate the individual cases. 


2. Collective modes - s wave ground state 


It can be seen from Eq. (1201) that in a pure angular 
momentum ground state (pure s or pure d), II LL = 0 if 
L ^ L' (due to orthogonality) and in non-complex order 
parameter ground state, II 12 = 0. Thus, in the s wave 
ground state, the only surviving bubbles are: 


n?i(Q) = v 2D 


= D 


—4L q + ( 1 — 


n 


( n 


V 2Ar 




-4L a - , 

9 ' 2Ar 


n&(n) 

n ff(n) = /o ->/ 2 (= Jo), ie (1,2) 


Jo(fi) 


(26) 


where L g = = ^ln-^r. This implies that (1) s and 

d channels are completely decoupled and (2) amplitude 
and phase sector are completely decoupled. As a result, 
the collective mode equation (det[Eq. (fl9l) = 0]) in the 
amplitude(phase) sector of angular momentum L reads 


TT LL 

ii ll(22) 


V L 


= 0, 


(27) 


where, V L is the interaction in the L th angular momen¬ 
tum channel (s or d). 
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Collective modes in the amplitude sector: The solu¬ 
tions to the mode equation are contained in 


1 - 


ft 

2Af 


J 0 (ft) = 0, 


1 - 


ft 

2Ar 


2 2 

/ 2 (ft) =-. 

Ud U s 


(28) 

(29) 


Neither of these equations has an undamped solution (see 
Appendix C). 

Collective modes in the phase sector: The solutions to 
the mode equation are contained in 



2Aq) W!i > = 0 ’ 

(30) 

2 2 

/ 2 (ft) = - . 

Ud U s 

(31) 


The first equation yields the soft BAG mode (ft = 0) 
and the second (arising from the d— wave sector) yields 
the well known BS mode d 2 ! 34 Looking Eq. a and Eq. 
m we see that the BS mode frequency is related to 
the competing T c values for the s and d channels. Thus 
detecting the BS mode at low temperatures and recording 
the T c of a sample gives direct quantitative estimate of 
the competing d- wave pairing interaction.— This feature 
will change for a multiband system. 

The above modes lie below the minimum gap in the 
system and hence are not damped. Notice that at bound¬ 
ary of the s and s + id state, where u s = Ud, the BS mode 
softens as expected. 


3. Collective Modes in the s+id ground state 


The collective modes in the s + id ground state are 
interesting. The surviving bubbles in this state are: 


(ft) = V2D 


(ft) = v 2D 


nff(ft) = V 2 D 


ng(ft) = V2D 


>/■ 


2A 


_ 2 /in _ + | a: _ 


ft 

2Ai 


2 A 

-2 / In— + al/ 2 (ft)- 
|A e | 


ft 

2An 


io(ft) 


/o (ft) 


2 / / “ 2 '"iii + (“--(4) 2 ) ,s(!i) 

■■f, 


-2 / fe In ^ + alh{Q) - 


ft 

2AJ 


nS (ft) = V2Da a adI 2 m, 


/a (ft) , 

(32) 


where a s ,d are to be found from T = 0 solutions of the 
gap equation as discussed in Appendix A. The non-zero 
nf 2 couples the s and d channels and also the ampli¬ 
tude and phase sector. This is expected from Eq. (l20l) 
because the ground state itself is a mixture of the two 
angular momentum channels and the order parameter is 


complex. The physical consequence of this non-zero bub¬ 
ble is that the 4 decoupled sectors now coalesce into two 
2x2 sectors formed out of the 

(1) s-phase and the d— amplitude components which 
yields 



ft 

2ZV 


a 2 s I 0 


V-d 1 2 


0(33) 


This sector contains the gauge-mode at ft = 0. 

(2) s-amplitude and the d— phase components which 
yields 





r 2 2 t 2 

I 0 - a s a d I 2 


= 0. 
(34) 


This contains the collective mode with mixed symme¬ 
try that adiabatically continues to the Bardasis-Schrieffer 
mode in the s-wave phase; we will refer to this mode 
henceforth as the mixed-symmetry Bardasis-Schrieffer 
mode (MSBS). This mixed symmetry mode is the analog 
of the d + id' and p + ip' clapping modes discussed in 
Refs. I 2 UI 23 H 24 I . 

The d— ground state can be similarly worked out (see 
Appendix C). There are no collective modes (other than 
the BAG mode) that propagate and hence we do not 
dwell on this further. 

Fig- Die) traces all the collective modes below the min¬ 
imum gap across the phase diagram. It is worth noting 
that all the collective modes could be found essentially 
from one mode equation. 


III. COLLECTIVE MODES IN A 3-POCKET 
MODEL 


We now move to the 3-pocket model which is more 
relevant for the FeSCs. Other than the multiband aspect, 
the approach is identical to the 1 band model. We shall 
thus focus on discussing the results and highlight the 
differences with the 1-band model. 


A. Model and phase diagram 


Here we study a prototypical FeSC system with one T 
centered hole pocket and two M centered electron pock¬ 
ets (the latter two are from the same band-and hence are 
related by symmetry). The model and the pairing inter¬ 
actions between the fermions is schematically shown in 
Fig. [2] The interactions in explicit form can be written 
as (only the leading harmonics are retained; ei —> +1 
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and e 2 —> — 1): 


v hh (k,p) = uz + ugf k f p , 

V eiei (k,P) = U s e +U d e ( 1 )( 1 ), 

V e2 e 2 {k,P) = U: + U d (- 1 )(- 1 ), 

y eie 2 (fc,p) = u s eie2 + u d eie2 { i)(-i), 

= £^e + £^ e / fc (l), 

v he2 (k,p) = U° he + U d Jk{-l), 

(35) 


where fk = \/2cos9k. The self consistency gap equations 
for this model read: 



A- = 

p 

A ei = 


A 62 = 



V^ r A(W(? + V^ ei A e . 1 W* 1 + (eiOe 2 ) , 

p,k k k pk k k v ' \ 

y ^? 1 a ( w ~1 , 

pk k k J 

v^Ajw^l, 

pk k k J 

(36) 


FIG. 2: The interactions in a 3 pocket (1 hole and 2 electron) 
model. The interactions in light grey font are unimportant as 
far as the main message of the work with applications to hole 
doped FeSC is concerned and are thus set to zero. 


where W~ = tanh, with x £ {h,e i,e 2 }. We fur- 

ther assume, for the sake of simplicity of presentation, 
that the electron and hole bands have identical disper¬ 
sions: resulting in the same density of states. This choice 
of interactions requires the gap structure to assume the 
form 


then need to define a few more dimensionless parameters 
in analogy with the 1-band case 


a 


X 

S 


r s 


AJ/Ao; a d = A^/A 0 , x e (ft, e), 

u i 'a — u i 


(39) 


A£ = A^ + A h d f p , 

Af = At + A e d , 

At* = A® — A e d . (37) 

To minimize the parameter space we set the following 
interactions to zero: U^, U t . The rationale behind 
this lies in the fact that we want to model a system driven 
to a d— wave state by the hole pockets (hence we retain 
U d ) and the d— wave character in the rest is induced due 
to the interband coupling terms (hence we retain U d e ). 
The reason behind the choice of this model is related to 
our desire to eventually address the Raman experiments 
on the hole doped FeSC, and is explained further in Sec. 
ED The s —wave character in this system driven by U^ e . 
For brevity, we introduce the following dimensionless con¬ 
stants for the interactions 


v s ,d = v 2 dU s ^, u d = u 2D U d . (38) 

Recognizing the v s and u d are the main ingredients for 
our problem and that v d is only needed to induce SC in 
the electron pockets, we set v d = zu d ; where the ratio 2 is 
set to some constant. Keeping v s fixed (as in the 1-band 
case), the system now has u d (and 2 if one so desires) 
as the tuning parameter of the model. This is sufficient 
to generate the s to s + id to d phase diagram. We will 


where Ao is the s-wave gap on the hole pocket when 
v d = 0. This model is easily solved at T = 0 and T = T C 
(see Appendix D). The results for the boundaries of the 
s + id state are given as follows: at T = T C1 the critical 
u d (< 0) is the solution to (larger \u d \ favors a d —wave 
state) 




V V s 


V2 



(40) 


There are two points marking the boundary of the s + 
id state at T = 0. The s-wave side boundary (setting 
A d ’ h —> 0 and requiring r d to be arbitrary) yields 


— - 2r s = 2v s In |r s |, 

r s 

crit.i “A, - y/rj +4z 2 
u d = - ^2 - 


(41) 


The d— wave side boundary (with A^ ,s —> 0, r s arbitrary) 
u™*’ 2 is the solution to 


1 

2v 2 s 

and r d satisfies 



r d 


— +11 


—— +c 2 — In \r d \ 

zu d 2 


zu d 


zu d (c 2 - In \r d \), 


(42) 


















where C 2 = J / 2 ln|/| = 0.153. These boundaries are 
shown in Fig. Ha). The detailed solution for the gap 
components as a function of u d presented in Fig. H^)- 


B. Collective modes 

For multiband systems, we follow the same procedure 
to derive the collective mode equation. Unless we con¬ 
sider Cooper pairing between different bands, there will 
be no interband nf/^s of the type / G h G e . Then, Eq. 
(119p will be generalized to matrix equation in the band 
space. Ilf/ is then evaluated in the relevant ground 
state. This generalization is given by: 

E {n& i[6a b ]-2[Vj; 1 U6 LL 'S i t}6F& = 0,(43) 

j,L',b 

where a, b are the band indices. The explicit form of 
this equation is discussed in Appendix E. The only off- 
diagonal elements (in the band space) arise from \V L \~ 1 . 
As before, we drop the coupling to the density channel 
and work with a 2 s / d x 2 real / imag x 2bands = 8x8 matrix 
space. The interaction matrices in each angular momen¬ 
tum channel are given by 

™=(° »)■ W1= (£ «)■ < 44 > 


1. Collective modes in the s-wave ground state 


As in the 1-band scenario, n sd and II 12 are zero due 
to the symmetry of the ground state. The non-zero ITs 
are given by 


n;;,/, = "2D 

nfl.e = V2D 

n S 22,h = V 2D 
= "2D 


- 4 L h + 1 - 


n 

2Ar 


7f(U) 


-8 L e + 2 (a®) 2 - 


n 

2Ar 




-4 L h ~ 

-8L h -2 


n 

2Ar 


O 

2a; 






(45) 


where L d e = b In ?he , ■ A dd involves changing I 0 I 2 , 

I I 

but in the s —wave ground state Iq = I-i- Also for x G 
(h, e), 


V2D ^O^m(f^n) = 


1 


(f*Y 


E* 


(E-) 2 +(fe) 


2 ’ 


{E X Y = U- +«) 2 + (^) 2 /f. (46) 


8 

For the s —wave ground state, sectors in the s and d 
channels decouple. Further, the amplitude and phase sec¬ 
tors decouple. This leads to the following four decoupled 
equations: 


n 

n 


n ii h 

l n 

n tth 

f n 


dd 
11, e 


dd 

22,e 


ss ttss 
ll,h Ai ll,. 


ss tT' ss 

22,h iA 22,. 


+ 


+ 


2 Ud 

v l 

2u d 


V 


,.2 ’ 


4 

772" ■ 


(47) 

(48) 

(49) 

(50) 


It helps to note that at T = 0 the self consistency equa¬ 
tions tell us that 

1 

2 L e = — -, 

2 r s v s 

2 L h = - r -E. (51) 

Vs 

Eq. (1471) is the s —amplitude sector which has no solution. 
Eq. gED is the s—phase sector. This sector contains the 
BAG mode and the damped Leggett mode (blue dots in 
Fig. Htop) where the real part of LHS = RHS). Using 
T = 0 relations, Eq. (|48|) gives 



4T/tdg + 4 L(,Iq + 



T h T e 
2 o 1 o 


0(-52) 


Since Iq > 0, the only solution is = 0 the BAG 
mode. The other solution, a Leggett resonance, has an 
imaginary part and is thus damped. Eq. (® is the 
d— amplitude sector which also has no solution. Eq. (EB 
is the d— phase sector which contains the BS mode. Note 
that there is one true mode (below the minimum gap) 
and a resonance in the continuum as shown by the blue 
dots in Fig. Hbottom). To see that this true generically 
in the model, we again substitute for Id’s and use T = 0 
relations to we find- 


where 



4 


T h T e 
1 q 1 o 



S = 
A = 


u d 

1 - r. - 


d 


V S V% 

1 

-b r 

2 r„ 


UdV s \ 




> 0 . 


(53) 


(54) 

(55) 


Comparing Eq. m and Eq. (1541) we see that in the 
s —wave phase S > 0. The mode equation can be cast 
into LQ, a + Afl 2 — S = 0. This clearly has two solutions 
with fl 2 > 0 and < 0. The U 2 > 0 solution is given by 
the implicit in U equation 


U 2 = 


-A + VA 2 + 4 LS 


A > 0, L > 0. 


(56) 


2 L 
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FIG. 3: (a) The phase diagram in the 3-pocket model, (b) The evolution of the gap components with Ud in the 3-pocket 
model at T = 0. The boundaries of the s + id phase are marked with arrows. Red(Blue) represents the hole (electron) pocket 
and the solid(dashed) line corresponds to s— wave(d—wave) component of the gap. The grey line at zero is the size of error in 
the numerical calculation due to the choice of grid and resolution parameters, (c): The (undamped) collective modes across 
the phase diagram in different channels. The dashed black line is 2A m i n in the system. The dots indicate the s + id phase 
boundaries at T = 0. In all the figures z = 1/2 and v s = 0.2. 


This is the only real solution, representing the well de- are: 
fined BS collective mode. As we approach the s + id 
boundary in the s— state, S —> 0+ and the mass of the 
collective mode approaches zero (as expected). Note 
at this stage the following differences with the 1 -band 
model: (1) We get a damped Leggett mode; (2) we get Ilff h 
two BS modes where one is damped; (3) The BS mode 
frequency is no longer related to parameters determined 
at T c because of temperature dependence of the gap ra- n n e 
tios. Depending on the magnitude of the effect (which 
depends on a chosen microscopic model), this should be 
important for quantifying the experiment with the model. ^ 22 ,h 


n 


22, e 


n 

n 


sd 
12 ,h 

sd 
12,e 


n dd 

11, h 


n 


dd 


I jdd 

ii 22,/i 


n 


dd 

22,e 


2. Collective modes in the s + id ground state 
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V2 D 
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—8L e + 2 ( a e s Y- 
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2Ar 


m 


-4 J^ Lh + (a h d ) 2 lZ(n)- l£(0) 


-8L e + 2(c^) 2 /f(D)-2 


n 

2Ar 


TO 


v 2D [TOTO)], 

V2D [2c^/ 2 e (0)], 


V2D 


V2D 


V2D 
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—8L e + 2 ( (a e s ) 2 — 


fl 

2Ar 


TO 


-4 j g f 2 L h + (a h d ) 2 I%(n) - l£(n) 


—8L e + 2(a^) 2 /|(fl) — 2 


Cl 
2A r 




( 57 ) 


As in the 1-band case, the s, d, amplitude and phase 
sectors get coupled due to Ilf 2 7 ^ 0 . The non-zero ITs 


where Lh^ e = 4 hi 


2A 


(As' e ) 2 +fil(Aj ,e ) 2 


The T = 0 solutions suggest the following useful rela- 
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Phase-s 




Phase-d 




0/A„ 


FIG. 4: Collective mode solutions in the multi-band s —wave 
ground state. Top: The s-wave order parameter phase fluctu¬ 
ation sector - solution to Eq. (1481) The blue dots indicate the 
frequency of the collective mode/resonance(if the imaginary 
part is non-zero). The f2 = 0 solution is the BAG mode; and 
the Q between the coherence peaks is the Leggett resonance 
(as in MgB; 2 ™“) . Bottom: The d- wave order parameter phase 
fluctuation sector - solution to Eq. (| 50 ll. The dots indicate 
the frequency of the collective mode/resonance. There is a 
conventional BS mode below the minimum gap and another 
resonance in between the two s—wave gaps. Here is used pa¬ 
rameters v s = 0.2, Ud = —0.18, z = 1/2. 

tions 

1 = -2v s r s [2L e ], 

1 “ -tJ™ 

1 = -u d J f 2 [2L h ] - 2v d r d [2L e \, 

1 = f f[2L h }. (58) 

Td J 

We thus get two decoupled 4x4 sectors formed out of 
22ss — 21 sd — 11 dd and llss — 12sd — 22 dd parts, of the 
form 


/ 

TTSS 

ii 22,/i 

_2_ 

v s 


0 \ 


_2_ 

v s 

TTSS 

Ai 22,e 

0 

TTSd 

AA 21,e 


yrsd 
AA 12 ,h 

0 

Tjdd 

AA 11 ,h 

_2_ 

Vd 

V 

0 

T\ sd 

iA 12,e 

_2_ 

Vd 

ridd 2 Ud 

iln > e + / 


FIG. 5: Solutions to Eq. (1591) (top) and its counterpart with 
Nambu components 1 and 2 exchanged(bottom) in the s + id 
phase. Solution to Eq. Eq. (1591) contains the BAG mode 
(blue dot at = 0) and the other equation has the mixed 
symmetry collective mode (blue dot at finite fl). This mode 
goes soft at the boundaries of the s + id state. Note that in 
the s + id state, phase and amplitude degrees of freedom are 
mixed in both equations. Here we used parameters v s = 0.2, 
u d = -0.231, z = 1/2. 

The other equation is given by (1 <+■ 2). The solutions 
are plotted in Fig. [5] We obtain the BAG mode and 
a single mixed symmetry collective mode. The result of 
tracking the collective modes (with no imaginary part) 
as function of u d is shown in Fig. [31(c). There are no 
other damped resonances. 


IV. DISCUSSION 

Having understood the spectrum of collective modes 
that can exist in a specific multiband model for a SC 
with competing s and d channels, we now discuss in 
some detail the practical motivation behind the choice of 
our model. As mentioned earlier, FeSC are complicated 
systems with a varying range of FS topologies and usu¬ 
ally the minimal model depends on the particular fam¬ 
ily of FeSC of interest. We chose to study the case of 
s to d transition as it is supposed to host the exotic 
s + id state whose detection is challenging. We pro¬ 
pose the identification of the MSBS collective mode to 
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detect the s + id state. Heavily hole or electron doped 
FeSC are the best candidates for the detection of the 
s + id state. Our model was designed keeping in mind 
Bai_ a; K a ;Fe 2 As 2 near x ~ 1. This is an interesting ma¬ 
terial which is known to be a fully gapped s± SC at 
x ~ 0.4^r—. It is also known to be nodal at x ~ 1. 
The symmetry of the latter compound is heavily de¬ 
bated: thermal conductivity^ data is interpreted in favor 
of d- wave which is in line with functional renormalization 
group prediction^. This is strongly contradicted by the 
photoemission dat a^— 2 which supports s— wave and also 
has theoretical support^. In either case, traversing from 
x ~ 0.4 to x ~ 1 may require the system to go through 
exotic states like s + id or s + is' as the case may be. The 
collective modes in the s + is' stat o 26 ’ 27 have Leggett-like 
collective modes (with s—symmetry) that soften at the 
transition. We have now shown that the collective modes 
in the s + id state have a mixed symmetry BS mode 
(with contributions in both s and d symmetry channels) 
which also softens at the boundary of the s + id state. 
The onset of the s + id state can thus be identified by 
simultaneous (and/or correlated) observation of the on¬ 
set of modes in both the s— and d— channels of Raman 
spectroscopy. Experimentally, there is strong indication 
of a mode in the d— channel as claimed in some recent 
works i 14 : 15 The signal in the s— channel is not conclusive 
which leaves us with two possibilities - (i) The system is 
still in the s-wave ground state, (ii) the screening effects 
in the s channel of Raman spectroscopy washes out the 
mode features, in which case the intensity of the s com¬ 
ponent in the MSBS must be carefully analysed. Both 
scenarios are unexplored theoretically as far as the ma¬ 
terial is concerned. 

Its worth mentioning that the proposed d— wave state 
in the x ~ 1 sample has the property that SC is driven 
by the hole pockets^. We thus chose to study the model 
where SC was driven by d— wave interaction within the 
hole pocket. 

We should warn the readers, however, that the exis¬ 
tence of collective modes is different from actually ob¬ 
serving them through experimental probes. One of the 
most promising probe to detect these modes (along with 
their symmetries) is electronic Raman spectroscopy. But 
the coupling of the modes to a Raman probe needs a spe¬ 
cial attention due to screening effects in the s— channel. 
The distribution of intensities across the s and d parts of 
the mixed symmetry collective modes is another aspect 
which requires careful attention. Thus, while a quanti¬ 
tative mapping of our results to Raman still needs more 
work, we are in a position to at least suggest the num¬ 
ber and nature of the collective modes to be expected in 
FeSC. 

Finally, although we know form the efforts of 
Ref. [l8| the possible modes in a ‘conventional’ nodal 
d— wave state, there are other d- wave models like fully 
gapped d— wave state proposed for systems like alkali 
intercalated-FeS o 44 ' 45 that need exploration in terms of 
the behavior of collective modes and is left for a future 


effort. It is also worth noting that as far as mixed sym¬ 
metry modes are concerned, such an analysis could be 
applied to systems under uniaxial strain where we gen¬ 
erate a similar s — d competition^ 


V. CONCLUSIONS 

To summarize, we have used a simple linear response 
approach to study collective modes in a SC with com¬ 
peting s and d channel instabilities. We work close to 
the region where s and d channels are nearly degenerate 
such that the system supports an exotic s + id state. We 
model a system that can be tuned from a pure s state 
to pure d state through the s + id state. Although our 
goal is the description of a multiband FeSC system, we 
first used a 1-band model to demonstrate the simplic¬ 
ity of the approach. In this 1 band model, we find the 
massless BAG mode throughout the phase diagram in 
the s— channel. In the d— channel, in the s— wave ground 
state we find the BS mode and in the s + id ground state 
find a mixed symmetry collective mode that softens at 
the boundaries of the s + id state. These boundaries 
were also calculated analytically. This mixed symmetry 
mode, which couples amplitude and phase sectors, oscil¬ 
lates with both s and d wave components and only exists 
in the T-broken s + id SC state. It is interesting to note 
that such a mode, which we refer to as a mixed-symmetry 
Bardasis-Schrieffer mode, corresponds, in a certain sense, 
to modes in T -breaking ground states that have been dis¬ 
cussed before, usually in situations where the two com¬ 
peting interactions correspond to degenerate basis func¬ 
tions of a 2D representation, e.g. “clapping modes” in 
p + ip' or d + id' situations'^ More generally if there 
are two competing representations, a mode of this type 
is possible^. 

We generalized our approach for a multiband system 
and found (1) the usual BAG mode in the s channel 
across the phase diagram; (2) damped Leggett mode be¬ 
tween the electron and hole gaps in the s— wave state; 
(3) 2 BS modes (where one is damped, like the Leggett 
mode) in the s— wave state; (4) A mixed symmetry col¬ 
lective mode in the s + id state that softens at the 
boundaries. Based on the multiband model that was 
designed to minimally reproduce the qualitative effects 
of Bai_ x Ka,Fe 2 As 2 near x ~ 1 , we suggest detection of 
the above ‘symmetry-selective’ collective modes. We pro¬ 
pose that such a systematic search can eventually settle 
the longstanding debate about the pairing symmetry for 
the x ~ 1 samples. 
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VI. APPENDIX 
A. Solving the 1-band model 


We solve the model described by Eq ([3]). For analytical answers will only look at T = 0 and T = T C points to study 
the special points in the phase diagram of the model. We conjecture (based on continuity) that there are no other 
special points in the phase diagram. At T = T C (A s ’ d —> 0), Eq. © leads to 

A s = —u s A s [2L c ], 

A d = -u d f /f A d [2L c ]. (60) 

•J 6 

where f g stands for the angular integral at the Fermi surface, L c = |ln|^. Using f g f 'i = 1 we see that 

2L C = min ; (61) 

l u s u d ) 

with the s + id state being realized when u d = u s . At T = 0, as discussed, there are three regions. For the s— wave 
state Eq. © yields: 


1 = —u s ln^, A d = 0. 

A s 

This requires u s < 0 (attractive) and also gives the T = 0 value for A s . We shall now define 

2A 1 

2 L = In—— =-. 

An Uc 


(62) 


(63) 


A§ is the s —wave gap in the model at T = 0 with no competing d —wave interaction. As the d —wave interaction ( u d ) 
grows, A d remains zero up to the point where the second equation in Eq. © first has a non-trivial solution, then A d 
begins to grow and A s begins to suppress. In the s + id state we have 


1 = —u s / In 

Jg 


2A 


A^ +A If I 


1 = ~u d /pin 


2A 


AH AH 


(64) 


In the notation of the main text, Eq. (1641) can then be cast into a simpler form- 
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-=2 L 0 


1 

u d 


I In— = 

6 yja% + ifa%f: 


-= 2 L n 




2^,2 f 2 
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« fe ^hh! 7 | 


(65) 


or, 
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u s 


lna s = 

1 

u d 


~i h V 1+??2/ l’ 

- f 0 {f l~ 1)ln \/ 1+rl2f l 


( 66 ) 


where 2 L g = ln^. This system of equations gives the gap ratios a s , d as a function of the parameter u d ■ 
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B. Explicit form for the 1-band mode equation 

In explicit form for the 1-band case the collective mode equation is 


/ Hff 

- 2 P - 1 

n i2 

TTSS 

aa 13 

TTSd 

AA 11 

TTSd 

aa 12 

n !s \ 


/ 

-6A r ’ s 


n|! 

n|| - 2 P - 1 

TTSS 

aa 23 

TTSd 

aa 21 

TTSd 

aa 22 

n 



SA 1 ’ 8 


TTSS 

Ai 31 

TTSS 

aa 32 

n II - V " -1 

s,q 

n If 

T\ sd 

aa 32 

TTSd 

li 33 



SD S 


T\ ds 

Ai ll 

T\ ds 

aa 12 

Tfds 

li 13 

nff - 217 1 

Tjdd 

aa 12 

vjdd 

aa 13 



-5A R ' d 


TJds 

aa 21 

T\ ds 

aa 22 

Tfds 

li 23 

Tjdd 

aa 21 

yrdd _ OT/ -1 
aa 22 z v d 

tt dd 
aa 23 



SA Rd 

V 

T\ ds 

aa 31 

ttcZs 

aa 32 

Tfds 

X1 33 

Y\dd 

aa 31 

rrdd 

aa 32 

ng - y 


l 

SD d 


Ignoring the density channel we get: 


( n?! 

- 2V- 1 

n i2 

TTSd 

AA 11 

N 


/ 

-SA R ’ S \ 


n 21 

nil - 2P- 1 

n|f 

nil 



SA Is 


TJ ds 

AA 11 

TJ ds 

aa 12 

n dd - 217 1 

tt dd 
aa 12 



-SA R ’ d 

V 

TJ ds 

aa 21 

TJ ds 

aa 22 

Tjdd 

AA 21 

n 22 - 2I7 1 J 


V 

5 A d ’ d / 


All information about collective modes are obtained by looking at the determinant of the above matrices. Depending 
on the relevant ground state, some of the IPs (as discussed in the main text) are zero, thereby simplifying the matrix 
structure. 


C. Finding solutions to the 1-band collective mode equation 


Here we present the graphical solutions to some of the equations presented in the main text. Fig. [ 6 ] shows the 
s —wave sector solutions (amplitude and phase) to the collective mode equation in the 1 -band s—wave ground state. 
Fig. [7]shows the d— wave sector solutions (amplitude and phase) to the collective mode equation in the 1-band s —wave 
ground state. Fig. [5]shows the collective mode solutions in the two coupled sectors in the 1-band s + id ground state. 

In the s + id ground state, the matrix equation in Eq. m becomes block diagonal in the ls- 2 d and 2 s-ld sectors. 
This mixes amplitude and phase and the s and d channels. The equations in the two sectors are: 


(n;f-A )( ng--)-<##£ = o, 

U s Ud 

(HI5- -)-aWA = 0, 

Us Ud 


which simplifies to 




Io - aWA 




2 2 r2 


0, 

0. 


(69) 


( 70 ) 


The final forms are given in Eqs. (l33l) and (13H) . 
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FIG. 6: (Left) Graphical solution to the amplitude sector Eq. (12811 in 1-band s —wave ground state with s —wave fluctuations. 
No solution for 11 < 2 Aq. (Right) Graphical solution to the phase sector Eq. (1301) with 1-band s —wave fluctuations. We get 
the BAG mode at 11 = 0. 







> 

- Re[LHS] 

- lm[LHS] 

-RHS 


■25-.-■-■-■-■- 

0 1 2 3 4 5 6 

^ /A 0 


FIG. 7: (Left) Graphical solution to the amplitude sector Eq. (1291) in 1-band s —wave ground state with d —wave fluctuations. 
No solution for SI < 2A§. (Right) Graphical solution to the phase sector Eq. (1311) with d —wave fluctuations. We get the BS 
mode at finite H that softens when Ud = u s - Here, u 3 = —1 and Ud = —0.6. 



^0 





FIG. 8: (Left) “Amplitude” sector Eq. (13311 - with H 2 removed in the s + id ground state. No solution. (Right) “Phase” 
sector Eq. (|34l) in the s + id state. There is only one solution indicated by the intersection of LHS and RHS. As Ud is tuned, 
this solution starts from 0 at the s —boundary and reaches 2A S at the other boundary-which in turn —y 0. 
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1. The d—wave ground state 


mi = V2D 


m s 2 = vw 


n if = v 2D 


II 22 — V 2D 


f-P 

[P 

Ik E 

f P 

li E 






— ) / 0 + a d I 2 


ft 


n , r 

— ) 1 2 + a d h 


The gap equation yields 


1 

u d 


Also, 


P , 2A 
2B = - ln ^ + C " 


Ik E u d 


These yield the mode equations: 




f2\ „ 

— ) / 0 - a d /2 


2 2 

—-1, 

Ud U s 

2 2 

Ud u s 


n 


h = 0 , 


(71) 


(72) 


(73) 


-y ) h~ a 2 dl 4 = 


It can be explicitly checked that, other than the BAG mode (12 = 0), there is no solution to these equations. 


(74) 


D. Solving the 3-pocket model 

Using the interaction form defined in Eq (15S1) and performing Jg — > V 2 D f y~ f de, the gap equations at T c can be 
written as 


A- = — 

p 


A 6 - 1 = - 

p 


A 6 - 2 = - 


Udfp J s f% + v d^Jp ^ A e J + Vs^Jp J s AA + ( ei o e 2 ) 
v.jA h + v d jf % A\ [2 L Tc ), 


[2 L Tc ], 


m, 

rn h 


v s [a h -v d I fa A 
Je Je 


[2£rJ, 


( 75 ) 
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where f deWS = 2Lt c = In In the pure s —wave state, we get: 


a? = 


AS = -u s 

1 


—Aj[2L To ], 

m R 


^ 2jLTc v'W AS 


In the pure d—wave state, we get: 


Ag = -u d A^[2L T J - 

AS = -u d 


[2A2][2Lr e ], 


=> 21/7’ = 


Ud 


—Ag[2L Te ], 

m K 


+ y/^d + 8^ Aj 


4u d 


4uj 


«d + -\AJ-i-8ufV 


The s—wave solution will be dominant when 


< 1 + 


/ Ud 


V2 V 1 


(76) 


(77) 


(78) 


with the equality reaching when the two solutions are degenerate. This will be the s + id point at T = T c . We shall 
now set m e = rrih to continue. Its going to be cumbersome otherwise. 

At T = 0, we get the following system of non-linear equations: 


where, 


As = 

-Us[2A'][2L e ], 

AS = 

-Us A $J[2L h ], 

and 


A d = 

-u d A h d f fpLy 

J 0 

II 

CD'S 

<1 

-v d A h d f fl[2L h ] 

J 0 


v d [2 AS] [2L e ], 


2 L h 

2 L e 


In 


In 


2A 

\J (A^ ) 2 + /|(A()) 2 ' 

2A 

\J (Af ) 2 + /J(AS ) 2 ' 


(79) 


(80) 


In the pure s— or d— phase we simply set the other gap to zero. At T = 0, there are two points: where Ad —> 0 (the 
s-wave boundary) and where A s —> 0 (the d —wave boundary). 


1. s—wave boundary 


At this boundary, we have 


A^ = -Us [2 AS] [2 L e ], 
AS = -us A h s [2L h ], 
and 

A h d (l + u d [2L h ]) = -v d [2 AS] [2L e ], 
AS = -v d A h d [2L h }. 


( 81 ) 
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The way to proceed is to realize that A d A —> 0 but A d /A d can be arbitrary. This arbitrariness can be removed by 
eliminating the A^’s from the last two equations. This will later give the constraint equation for u d to get a non-trivial 
d —wave component. The first two equations should be used to compute the A h J e . These yield 

2 


-t u d 

1 ^ r. —, 


where r s = 


and r s satisfies 

—— 2 r s 
r„ 


K 


= 2v a In |r s |. 


(82) 


Note that at T = T r , r. 


= -75- AtT = 


this will obviously depend on v s . Recall that v d = zu d and u d is our 


tuning parameter for the phase diagram. This yields the critical u d to be 


cr it _ 
u d ~ 


~ r s ~ \fif 

2z 2 


4 z 2 


(83) 


And r s is to be solved as before. The result is presented in Fig. [31(a) . 


2. d—wave boundary 


At this boundary, we have 


A h s = -v a [2At][2L e ], 

A e s = -v s A h s [2L h -c 0 ], 
and 


A h d {l + u d [2L h -c 2 ]) = -v d [2A e d ][2L e ], 

A e d = -v d A h d [2L h -c 2 ], (84) 

where, Co = / In |/| = —0.347, and c 2 = J f 2 In |/| = 0.153, such that c 2 — cq = \. We then follow the same logic as 
for the s—wave boundary and eliminate A) 6 and get 


= [2L h -c 2 + ^][2L h -]n\r d \], 


of r d 

2Lh — c 2 =-, 

(85) 

zu d 

r d satisfies 


-2 r d = zu d (c 2 — In |fd|), 

(86) 


where z = — and r d 

'll. A 


-rjt. We just find the u d for which Eq. [85] is satisfied. The result is presented in Fig. ETa). 

“d 


E. The explicit form of the 3-pocket mode equation 


/ 

T-fSS 

Ai ll ,h 

0 

ttss 

AA 12 ,h 

0 

TTSS 

aa 13 ,h 

TT sd 

li ll ,h 

0 

Tfsd 

il 12,/i 

0 

TTsd 

AA 13,/i 


0 

yrss 

AA ll,e 

0 

TTSS 

AA 12,e 

TTSS 

AA 13,e 

0 

TT sd 

li ll,e 

0 

fTSd 

ii 12,e 

TTSd 

AA 13,e 


T-TSS 

iA 21 ,h 

0 

TTSS 

aa 22 ,h 

0 

TTSS 

aa 23 ,h 

TT sd 

1J -2l ,h 

0 

Tfsd 

X1 22 ,h 

0 

TTsd 

AA 23,h 


0 

TTSS 

AA 21,e 

0 

TTSS 

AA 22,e 

nr 3 ,e 

0 

TTSd 

11 21,e 

0 

ii 22,e 

niL 


T-TSS 

AA 31,h 

TTSS 

AA 31,e 

TTSS 

aa 32 ,h 

TTSS 

AA 32,e 


TT sd 

lx 3l ,h 

TT sd 

li 31,e 

Usd 

il 32,/i 

fTsd 

il 32,e 

Mg 


T\ ds 
AA 11 ,h 

0 

T\ ds 
aa 12 ,h 

0 

TTf/s 
aa 13 ,h 

Y\dd 

LL ll,h 

0 

rjdd 
AA 12 ,h 

0 

T\sd 

1 A 13,/t 


0 

TTC^S 

AA ll,e 

0 

T\ ds 

AA 12,e 

T\ ds 

AA 13,e 

0 

it d d 

AA ll,e 

0 

Y\dd 

AA 12,e 

TTSd 

AA 13,e 


Y\ ds 
aa 21 ,h 

0 

T\ ds 

aa 22 ,h 

0 

TTds 

aa 23 ,h 

rjdd 

LL 21,h 

0 

tt dd 
aa 22 ,h 

0 

tt dd 
aa 23 ,h 


0 

TTCte 

AA 21,e 

0 

T\ ds 

AA 22,e 

n 

0 

rjdd 

AA 21,e 

0 

T~[dd 

AA 22,e 

rjdd 

AA 23,e 

V 

T\ ds 

aa 31,Ji 

TTds 

AA 31,e 

T\ ds 

aa 32 ,h 

TTSS 

AA 32,e 

M d i 

rjdd 

LL 31,h 

tt d d 

AA 31,e 

yrdd 
aa 32 ,h 

rjdd 

AA 32,e 

M dd 


/ -6*2" 5 

-SA^ S 

K ' 8 

5Al’ s 
SD S 
-5A g d 
-8A*' d 
SA 1 / 
6A 1 / 
SD d 


( 87 ) 
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where M$ 3 L = 



/ 0 


0 

0 

0 

0 

0 

0 

0 

0 \ 

-1 

/ 

sa£’ s 


v. 

0 

0 

0 

0 

0 

0 

0 

0 

0 



-6A*' S 


0 

0 

0 

v s 

0 

0 

0 

0 

0 

0 



K ’ 8 


0 

0 

v a 

0 

0 

0 

0 

0 

0 

0 



6Al’ s 


0 

0 

0 

0 

2 vs 

0 

0 

0 

0 

0 



6D S 

2 

0 

0 

0 

0 

Q 

0 

u d 

v d 

0 

0 

0 



-5A*' d 


0 

0 

0 

0 

0 

v d 

0 

0 

0 

0 



-8Ay 


0 

0 

0 

0 

0 

0 

0 

U d 

V d 

0 



8A I h ’ d 


0 

0 

0 

0 

0 

0 

0 

v d 

0 

0 



h 

5A I ’ d 


V o 

0 

0 

0 

0 

0 

0 

0 

0 

2 Vi ) 

V 

e 

SD d 


( 88 ) 







